Experimental Schmidt Decomposition and State Independent Entanglement Detection 
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We introduce an experimental procedure for the detection of quantum entanglement of an un- 
known quantum state with a small number of measurements. The method requires neither a priori 
knowledge of the state nor a shared reference frame between the observers and can thus be regarded 
as a perfectly state independent entanglement witness. The scheme starts with local measurements, 
possibly supplemented with suitable filtering, which essentially establishes the Schmidt decomposi- 
tion for pure states. Alternatively we develop a decision tree which reveals entanglement within few 
steps. These methods are illustrated and verified experimentally for various entangled states of two 
and three qubits. 



PACS numbers: 03.67.Mn 

Introduction. — Entanglement is the distinguishing fea- 
ture of quantum mechanics and it is the most important 
resource for quantum information processing [TJ [2] . For 
any experiment it is thus of utmost importance to easily 
reveal entanglement, best with as little effort as possible. 
Common methods suffer from disadvantages. On the one 
hand, employing the Peres-Horodecki criterion [3] [4] or 
evaluating entanglement measures, one can identify en- 
tanglement in arbitrary states, however, it requires full 
state tomography. On the other hand, various entan- 
glement witnesses [414TQ] can be determined with much 
fewer measurements but give conclusive answers only if 
the state under investigation is close to the witness-state, 
i.e., they require a priori knowledge. 

Recently, it has been shown that the existence of en- 
tanglement can be inferred from analyzing correlations 
between the measurement results on the subsystems of a 
quantum state. Only if the state is entangled, the prop- 
erly weighted sum of correlations will overcome charac- 
teristic thresholds [11]. Here we further develop this ap- 
proach to obtain a simple and practical method to detect 
entanglement of all pure states and some mixed states by 
measuring only a small number of correlations. Since the 
method is adaptive it does not require a priori knowl- 
edge of the state nor a shared reference frame between 
the possibly remote observers and thus greatly simpli- 
fies the practical application. We describe two schemes 
where the first one essentially can be seen as a direct im- 
plementation of Schmidt decomposition, which identifies 
the maximal correlation directly. For bipartite pure sys- 
tems, this approach can conceptually be divided into two 
stages: (i) calibration that establishes the experimental 
Schmidt decomposition [12] [13] of a pure state by local 
measurements and suitable filtering and (ii) two correla- 
tion measurements to verify the entanglement criterion. 
The second scheme shows how to use a decision tree to 
obtain a rapid violation of the threshold thereby identi- 
fying entanglement. 

Entanglement criterion. — For a two-qubit quantum 



state p, Alice and Bob observe correlations between their 
local Pauli measurements 0& and 07, respectively. They 
are defined as the expectation values of the product of 
the two measurements, Tki = Tr[p(ak (S> 07)], with the so 
called correlation tensor elements Tki G [—1,1]. The lo- 
cal values Tfco (To/), with O being the identity operator, 
form the local Bloch vector of Alice (Bob). Using these 
measurements a sufficient condition for entanglement can 
be formulated as (TTJ [14] : 

^2 T ki > 1 => P is entangled. (1) 

k,l=x,y,z 

For pure states this is also a necessary condition, while 
for mixed states care has to be taken. For those the 
likelihood for detecting the entanglement decreases with 
purity. An extension of ([I]) can generally identify en- 
tanglement of an arbitrary mixed state, however, then 
loosing the state independence [11] . Note two important 
facts. First, Eq. can be seen as a state independent 
entanglement witness derived without any specific family 
of entangled states in mind. Second, to test whether the 
state is entangled, it is sufficient to break the threshold, 
i.e., it is neither required to measure all correlations nor 
to compute the density matrix [25]. Rather, it is now the 
goal to find strategies which minimize the number of cor- 
relation measurements. We show how this can be done 
by a particularly designed decision tree, or by identifying 
a Schmidt decomposition from local results and filtering 
when necessary. 

Schmidt decomposition. — Consider pure two qubit 
states. Any such state has a Schmidt decomposition 

tys) = cos 0\a) \b) + sintf \a ± ) \b ± ) , G [0, \ ], (2) 

where the coefficients are real and the local bases 
{\a} , \a±)} and {\b} , \b±)} are called the Schmidt bases. 
Once the bases are known, Alice constructs her local mea- 
surements g z i = I a) (a\ — \a±) (a±\ and a y > = i \a±) (a\ — 
i I a) (aj_|, and so does Bob in analogy. They can now 
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detect entanglement with only two correlation measure- 
ments because T 2 , + T 2 , , = 1 + sin 2 20 > 1 for all pure 

z z y y 1 

entangled states. Note, the laboratories are not required 
to share a common reference frame. 

In order to extract the Schmidt bases from experimen- 
tal data one starts with local measurements, determining 
the local Bloch vectors a (/?) of Alice (Bob) (those ele- 
ments are rela ted to the corre lation tensor coefficients 
via oti = T iQ I yj T 2 + T 2 + T 2 ) . We consider two cases. 

First, suppose that a pure state has non- vanishing local 
Bloch vectors. Their directions define the Schmidt bases 
of Alice and Bob up to a global phase <p. Writing these 
bases in the computational basis 



\a) 
\a±) 

\b) 
\b±) 



cosU\0) + e iVA sin U | 1) , 
sin^O)-^ cos^|l), 
cos£ B |0) + e^ B sin£ B |l), 
e^(sin^ s |0)-e^cos^ B |l)), 



(3) 



we see that the required coefficients can be in- 
ferred directly from the local Bloch vectors, a = 
(sin 2^a cos (fA, sin 2£a sin (£>a-> cos 2£a) on Alice's side, 
and similarly for Bob. The global phase of|&_i_) shows up 
as the relative phase in the decomposition i.e. \ips) 

cos (9 | a) 1 6) + sinfle^ \a±) \b±) (with \b±) = 



can be determined, for example, from the T yy 



as cose 
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b±)). It 
correlation 
If Bob would use 
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the basis {\b) , to build his observables o z n and a, 
the corresponding correlations T y i y n = sin 2d cos <p would 
vanish for cos(/> = and the two measurements T z i z n 
and Ty/yff would not suffice to detect entanglement. In 
such a case, however, the other two correlations, T x > y n 
and T y i x n , are non-zero, and can be used to reveal entan- 
glement. Therefore, the determination of <fi in the cali- 
bration is not essential if one accepts possibly one more 
correlation measurement. 

Second, in case of vanishing local Bloch vectors, the 
pure state under consideration |-0 m ) is maximally entan- 
gled and admits infinitely many Schmidt decompositions. 
In order to truly prove entanglement Bob can thus freely 
choose some basis, say computational basis, for which the 
state will now be of the form \i/j m ) = ~^(\ a ) |0) + |aj_) |1)). 
The basis of Alice can be found after filtering by Bob in 
his Schmidt basis: F = |0) (0| + e |1) (1| (for an actual 
implementation see experimental section). When Bob 
informs Alice that his detector behind the filter clicked, 
the initial state becomes 



(1®F)|VU 



= (|a)|0)+ e |oL)|l)). (4) 



Note that due to filtering a nonvanishing local Bloch vec- 
tor emerges for Alice. Thus the respective Schmidt basis 
can be found with the method described above and be 
used for the evaluation of T 2 , r + T 2 , 7/ . 

Decision tree. — Our second algorithm for entangle- 
ment detection does not even require any calibration and 
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Figure 1: The decision tree for efficient two-qubit entangle- 
ment detection. No shared reference frame is required be- 
tween Alice and Bob, i.e. they choose their local x,y,z direc- 
tions randomly and independently, which effectively gives rise 
to a basis {xa,Va, za} for Alice and {xb,Vb, zb} for Bob (not 
detailed in the Figure or the main text). The scheme starts 
with measuring T zz and follows at each step along the dashed 
arrow if the modulus of correlation is less than \ and along the 
continuous arrow otherwise. The algorithm succeeds as soon 
as ^ Tfj > 1 . The measurements in blue shaded area suffice 
to detect all maximally entangled pure states with Schmidt 
bases vectors along x,y or z. 



also applies directly to mixed states. Alice and Bob 
choose three orthogonal local directions x,y and z in- 
dependently from each other and agree to only mea- 
sure correlations along these directions. In Fig. [T] we 
show exemplarily which correlations should be measured 
in order to detect entanglement in a small number of 
steps. Starting with a measurement of T ZZ1 one contin- 
ues along the solid (dotted) arrow, if the correlation is 
higher (lower) than some threshold value (e.g. 1/2 in 
Fig. [I]). The tree is based on the principle of correlation 
complementarity [T5HT8] : in quantum mechanics there 
exist trade-offs for the knowledge of dichotomic observ- 
ables with corresponding anti-commuting operators. For 
this reason, if the correlation \T ZZ | is big, correlations 
\Tzx\-, \T zy \, \T XZ \ and \T yz \ have to be small because their 
corresponding operators anti-commute with the operator 
o~ z ®°~z- Therefore, the next significant correlations have 
to lie in the xy plane of the correlation tensor and thus 
the tree continues with a measurement of the T yy cor- 
relation. This concept can be generalized to multiqubit 
states, a decision tree for three qubits is given in the Ap- 
pendix. The number of detected states grows with the 
number of steps through the decision tree. Since con- 
dition (1) is similar to the purity of a state, obviously, 
the scheme succeeds the faster the more pure a state is 
(see Appendix for detailed analysis). Varying the thresh- 
old value does not lead to any significant changes in the 
statistic of detected states. 

Finally, we connect both methods discussed here for 
the analysis of multiqubit states. A numerical simula- 
tion for pure states reveals that the correlation measure- 
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Figure 2: Scheme of the experimental setup. The state 
|^r) = ^(|if) \H) + e iS \V) \V)) is created by type I SPDC 
process. An Yttrium Vanadate crystal (YVO4) is used to ma- 
nipulate the phase 5 of the prepared state. For preparation 
and analysis of the state half- (HWP) and quarter (QWP) 
waveplates are employed. Brewster plates (BP) can be intro- 
duced to make the state asymmetric and to perform the filter 
operation, respectively. 



merit along local Bloch vectors gives correlations close to 
the maximal correlations in more than 80% of the cases. 
Therefore, these local directions give an excellent starting 
point for the decision tree. 

Experiment. — For the demonstration of these new sim- 
ple analysis methods we first use two photon polarization 
entangled states. In the following we will thus replace 
the computational basis states by horizontal (|0) — » |if)) 
and vertical (|1) — > \V)) linear polarization, respec- 
tively. The photon source (see Fig. 2| is based on 
the process of spontaneous parametric down conversion 
(SPDC) using a pair of crossed type I cut /3-Barium- 
Borate (BBO) crystals pumped by a CW laser diode 
at a wavelength of \ pU mp — 402nm with linear polar- 
ization of 45°. It emits pairs of horizontally and ver- 
tically polarized photons which superpose to the state 
I*) = 4= (| if) |if) + e i5 \V) \V)) [19]. The spectral band- 
width of the photons is reduced to 5nm using interference 
filters and two spatial emission modes are selected by 
coupling the photon pairs into two separate single mode 
fibers. 

For the purpose of preparing any pure two-qubit state, 
the polarization of each photon can be rotated individu- 
ally by a set of quarter- (QWP) and half (HWP) wave- 
plates in each mode. By tilting an Yttrium Vanadate 
crystal (YVO4) in front of the BBOs, the relative phase 
S between the photon pairs can be set. Additionally, 
the state can be made asymmetric by removing a por- 
tion of vertically polarized light in one spatial mode with 
a Brewster plate (BP). In the last step of the experi- 
ment, the polarization of each photon is analyzed with 
additional half- and quarter waveplates and projection 
on I if) and \V) using a polarizing beam splitter (PBS). 
The local filtering of a maximally entangled state can be 
accomplished by placing a Brewster plate in front of the 
analysis waveplates. This Brewster plate reflects with a 
certain probability vertically polarized photons and to- 
gether with detection of a photon behind the Brewster 
plate implements the filtering operation Q. Finally, the 
photons are detected by fiber-coupled single photon de- 
tectors connected to a coincidence logic. 



Experimental Schmidt decomposition. — Let us con- 
sider the state shown in figure [3^1). The protocol starts 
with Alice and Bob locally measuring the polarization of 
the photons enabling them to individually determine the 
local Bloch vectors. For high efficiencies, e.g., possible in 
experiments with atoms or ions, the local measurements 
can indeed be done independently [26 j. If nonvanishing 
local Bloch vectors can be identified, one can proceed to 
the next step. For the example here, the local expecta- 
tion values are close to zero and filtering has to be ap- 
plied. By using a Brewster plate in front of Bob's analysis 
waveplate and under the condition that the filtering op- 
eration is successful, local Bloch vectors emerge (see Fig. 
[3]o) [27 , in this case we obtain T i = (0.000, 0.040, 0.334) 
and T k0 = (0.188,-0.034,0.336). 

In the next step Alice and Bob use their local Bloch 
vectors to realign their analyzers to the new local 
Schmidt bases {|a), \a±)} and {|6), |6j_)}, respectively. 
This diagonalizes the correlation tensor as depicted in 
figure [3^). Therefore, it is only necessary to measure 
T z > z n = 0.922±0.015 and T y , y n = -0.864±0.015 to prove 

entanglement since T?, „ + T?.„ = 1.597 ± 0.038 > 1. 

z z y y 

Hence, 2x3 local measurements are needed in the first 
step of the algorithm, if necessary three combined mea- 
surements are needed for filtering, and finally only two 
correlation measurements have to be performed for en- 
tanglement detection. 

Application of the decision tree. — In order to demon- 
strate the application of the decision tree we will apply it 
to three states. For the first state 4^(1 if ) |if) + \V) \ V)), 
whose correlation tensor is depicted in |4^i) the deci- 
sion tree (see Fig. [I]) starts with the measurement of 
the correlation T zz — 0.980 ± 0.015 and continues with 
T yy = —0.949 ±0.015. These two measurements already 
prove entanglement since T% z + Tjj y = 1.869 ± 0.041 > 1. 
For a second state, ^=(\R) \R) -\-i\L) \L)), (see Fig. [Z]d), 
we obtain a correlation of T zz = —0.056 ±0.015, close to 
zero. Consequently, the next steps according to our al- 
gorithm (Fig. IT]) are to determine the correlation T yy = 
0.978 ± 0.015 followed by T xz = -0.959 ± 0.015 , with 
their squares adding up to a value of 1.879 ± 0.041 > 1 
and hence proving entanglement. As a last example we 
consider the initial state of Fig. [3j According to our 
decision tree we need to measure T zz = 0.768 ±0.015, 
Tyy = 0.018 ± 0.015 and T yx = -0.922 ± 0.015, thus giv- 
ing a value of 1.440 ±0.036 > 1 and proving entanglement 
with only three steps. 

Many qubits. — For the demonstration of multiqubit en- 
tanglement detection, we use two three-photon polar- 
ization entangled states: the W state [20] and the G 
state [21 (see Fig. |5|. In order to observe these states, 

a collinear type II SPDC source together with a linear 

(2) 

setup to prepare the four-photon Dicke state D\ is used 
[22j [23]. Once the first photon is measured to be ver- 
tically polarized, the other three photons are projected 
into the W state. Similarly, the three-photon G state 
is obtained if the first photon is measured to be ±45° 
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Figure 4: Correlation tensors and density matrices of the 
experimental realization of two different states. The imagi- 
nary parts of the density matrices are negligible and therefore 
skipped. Using the decision tree, only the blue shaded corre- 
lations have to be measured for detecting entanglement. The 
errors of the correlations are < 0.015 for a) and < 0.023 for 

b). 



Figure 3: Demonstration of Schmidt decomposition of an 
maximally entangled state prepared in unknown bases. The 
correlation tensor and corresponding density matrix are de- 
picted for the unknown state a), after applying local filtering 
b) and for the state analyzed in the Schmidt bases c). It is 
important to note that only the blue shaded elements of the 
correlation tensors will be measured as this suffices to prove 
entanglement. The full correlations and the density matrices 
of the corresponding states are only shown for completeness 
and didactical reasons. 



polarized. 

The protocol for entanglement detection starts with 
observers locally measuring the polarization of the pho- 
tons enabling them to individually determine the lo- 
cal Bloch vectors. For the G state we obtain T^oo = 
(0.636,-0.008,-0.015), T 0j0 = (0.623,-0.092,0.010) 
and Took = (0.636,0.070,0.022). The local Bloch 
vectors suggest that the correlation T xxx is big. 
Therefore the decision tree starts with the measure- 
ment of T xxx = 0.904 ± 0.025 and continues with 
T xzz = -0.578 ± 0.025 (see Fig. 7 in the Ap- 
pendix). These two measurements already prove en- 
tanglement because T xxx + T xzz = 1.152 ± 0.038 > 
1. For the W state, the local Bloch vectors: T^oo = 
(0.016,-0.070,0.318), T 0j0 = (-0.010,-0.073,0.308) 
and Took = (-0.011,-0.0547,0.319) suggest that now 
the correlation T zzz is big. Indeed, we observe T zzz = 
—0.882 ± 0.025. The decision tree is the same as above 
but with local axes renamed as follows x — >• z — >• y — >• x. 
Therefore, the second measurement has to be T zyy . With 
T zyy = 0.571 ± 0.025 we again prove entanglement as 
T* xx +T 2 zyy = 1-104 ± 0.037 >1. 

Conclusions. — We discussed and experimentally im- 
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Figure 5: Density matrices of the experimental realization of 
the G and W state. The corresponding fidelities are equal to 
92.23% and 89.84%. 



plemented two methods for fast entanglement detec- 
tion for states about which we have no a priori knowl- 
edge. They are well suited for quantum communication 
schemes as the parties do not have to share a common ref- 
erence frame, making the scheme insensitive to a rotation 
of the qubits during their transmission to the distant lab- 
oratories. The two methods use a particularly simple and 
practical entanglement identifier [11]. One of them can 
be seen as experimental Schmidt decomposition and the 
other establishes a sequence of correlation measurements 
leading to entanglement detection in a small number of 
steps. 
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dation for Polish Science. CS thanks QCCC of the Elite 
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Appendix: Decision Tree 

This section is devoted to study the efficiency of the 
decision tree algorithm described in the main text. Some 
results in this section are analytical and some are nu- 
merical. In all our numerical investigations (unless ex- 
plicitly stated otherwise) we used the decision tree of 
the main text (for two qubits) and in cases when going 
through the whole tree did not reveal entanglement we 
augmented it with additional measurements of those cor- 
relations which were not performed until that moment. 
The order of the additional measurements also results 
from the correlation complementarity (anti-commutation 
relations) [15 . With every remaining measurement we 
associate the "priority" parameter 

= ^Pi j (T kj ) + Y / Pi j (T u ), (A.l) 

that depends on the measured correlation tensor elements 
of the decision tree in the following way 

( v _ J T^ n if T mn was performed before, 
^Wmnj-| o else. 

(A.2) 

According to the correlation complementarity there is a 
bigger chance that this correlation is significant if the 
value of the corresponding parameter is small. Therefore, 
the correlations with lower values of are measured 
first. 

Let us consider the following example. The measured 
correlations of the decision tree are as follows: T zz = 0.7, 
T xx = 0.1, and T yy = 0.4. Therefore, P xy = P yx = 
?L + T 2 yy = 0.17 , P xz = P zx = Tl x + Tl = 0.5, and 
P zy = P zy = T% z + Ty y = 0.65. Accordingly the order of 
the remaining measurements is as follows: first measure 
xy, then yx, next xz, zx, yz and zy. 

1. Two qubits 



This is because only three elements contribute to the cri- 
terion. Note that the value of the decision parameter 
is irrelevant here. Furthermore, a random choice of lo- 
cal coordinate directions does not help. Although more 
steps would be involved in the decision tree the sum of 
squared correlations is invariant under local unitary op- 
erations and therefore only for p > entanglement is 
detected for the Werner state. 



b. Entangled state mixed with colored noise 

An exemplary class of density operators for which the 
decision tree detects all entangled states is provided by: 

7=p |V;->(^-| + (l-p)|01)(01|, (A.4) 

where entanglement is mixed with colored noise 1 01) 
bringing anti-correlations along local z axes. This state 
has the following non-vanishing elements of its correla- 
tion tensor T xx = T yy = —p and T zz = — 1. Therefore, 
the decision tree allows detection of entanglement for this 
class of states in two steps. Note that the state is entan- 
gled already for an infinitesimal admixture of the Bell 
singlet state as can be shown for a random choice of lo- 
cal coordinate systems. 



c. Random mixed states 

Fig. [6] shows how the efficiency of the algorithm grows 
with the purity of tested states. The efficiency is mea- 
sured by the fraction of detected entangled states ob- 
tained from extensive Monte Carlo sampling in the two 
qubit state space. For nine steps the algorithm detects 
all the pure states. This is expected because Eq. (1) of 
the main text is a necessary and sufficient condition for 
the detection of entanglement of pure states. 



a. Werner states 

As an illustration of how the decision tree works for a 
well-known class of mixed states we first consider Werner 
states. It turns out that not all entangled states of the 
family can be detected. 

Consider a family of states 



p = p\ip )(ip | + (1 



■P)\l, 



(A.3) 



where \^~) = ^j(|01) - |10)) is the Bell singlet state, 

|l describes the completely mixed state of white noise, 
and p is a probability. Its correlation tensor, written in 
the same coordinate system for Alice and Bob, is diagonal 
with entries T xx = T yy = T zz = —p, arising f rom the con- 
tribution of the entangled state. The states (A.3) are en- 
tangled if and only if p > | , whereas the decision tree re- 
veals that these states are entangled for p > ^= w 0.577. 



2. Many qubits 

a. Algorithm for generation of the tree 

The principle behind the decision tree is the correla- 
tion complementarity [? ]. Correlation complementar- 
ity states that for a set of dichotomic anti-commuting 
operators {ai, . . . , a^} the following trade-off relation is 
satisfied by all physical states: 



T < 1 



(A.5) 



where T ai is the average value of observable ol\ and so 
on. Therefore, if one of the average values is maximal, 
±1, the other anti-commuting observables have vanishing 
averages. This motivates taking only sets of commuting 
operators as different branches of the decision tree. 
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Figure 6: Efficiency of the decision tree for two qubit random 
mixed states. 



b. Decision tree for three qubits 

This tree is an example of the application of the al- 
gorithm presented in the previous section. Fig. [7] shows 
only one branch of the whole tree. 

A numerical simulation reveals that the correlation 
measurement along local Bloch vectors gives correlations 



close to the maximal correlations of a pure multi-qubit 
state in more than 80% of the cases. Therefore, these 
local directions give an excellent starting point of the de- 
cision tree. 

Exemplarily, the branch begins with T xxx assuming 
that this correlation is big. If the local Bloch vectors 
indicate correlation along a different directions it is ad- 
vantageous to correspondingly change the elements of the 
decision tree (see main text). 
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Figure 7: One branch of the decision tree for three qubits. 
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